set more off

local J `1'
local sigma `2'
local sigmax `3'


forvalues i = 1/`J' {
gen x_`i' = `sigmax'*rnormal()
gen z_`i' = `sigma'*rnormal()
gen epsilon_`i' = -ln(-ln(runiform())) 
}
